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Abstract: We extend our earlier study of spin correlations in the decay chain D — > 
Cq, C — > Bl nea,r , B — > Al iar , where A,B,C,D are new particles with known masses but 
undetermined spins, Z near and £ far are opposite-sign same- flavour charged leptons and A is 
invisible. Instead of looking at the observable 2- and 3-particle invariant mass distributions 
separately, we compare the full three-dimensional phase space distributions for all possible 
spin assignments of the new particles, and show that this enhances their distinguishability 
using a quantitative measure known as the Kullback-Leibler distance. 
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1. Introduction 

In the recent paper [1], to which we refer the reader for motivation, notation and relevant 
references, we examined the distinguishability of different spin assignments in the decay 
chain D -» Cq, C -» Bl ncax , B -► Al fav , where A,B,C,D are new particles with known 
masses but undetermined spins, P ear and l far are opposite-sign same-flavour charged leptons 
and A is invisible. This was done by comparing separately the invariant mass distributions 
of the three observable two-body combinations: dileptons (run), quark- or antiquark-jet 
plus positive lepton (rriji+), and jet plus negative lepton (m^-). 1 

If P(m\S) represents the normalized probability distribution of any one of these three 
invariant masses predicted by spin assignment S, and T is the true spin configuration, then 
a measure of the improbability of S is provided by the Kullback-Leibler distance 



In particular, the number N of events required to disfavour hypothesis S by a factor of 
1/R under ideal conditions, assuming equal prior probabilities of S and T, would be 



By ideal conditions we mean isolation of the decay chain with no background and perfect 
resolution. Therefore N sets a lower limit on the number of events that would be needed 
in real life. The results for R = 1000 are shown in tables 1-3, reproduced for convenience 
from [1], where a discussion of them can be found. Recall that the notation used is DCBA 
with F for fermion, S for scalar, V for vector, so that squark decay in SUSY is SFSF and 
excited quark decay in UED is FVFV. Mass spectra I and II are SUSY- and UED-like 
respectively (see [1] for details). 

lr The three-body invariant mass rrijii was also studied but this is not independent of the two-body masses. 




(1.1) 
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2. Three-dimensional analysis 



To extract the most information from the data we should compare the predictions of 
different spin assigments with the full probability distribution in the three-dimensional 
space of mu, m ji+ and m ji-- The ambiguity between near and far leptons means that this 
given by 

P '(mu, m jl+ , mjl -) = ^f q [P 2 (m u ,m jl +,m jl -) + P 1 (mu,m jl -,m jl +)\ 

+ \fq [ p i( m ih m ji+i m ji-) + p 2(mu,m jl -,m jl +)\ , (2.1) 

where f q and f q = 1 — f q are the fractions of quark- and antiquark-like objects D initiating 
the decay chain and we use Pi t 2(mu,m'jf a ' r ,m i jf) on the right-hand side, assuming both 
leptons are left-handed, otherwise f q and f q are interchanged. The subscripts 1 and 2 refer 
to processes 1 and 2 defined in [1] and the factors of one-half enter because P\ :2 are both 
normalized to unity. 

Instead of trying to evaluate the three-dimensional generalization of the integral in 
eq. (1.1) analytically, it is convenient to perform a Monte Carlo integration. If we generate 
mu, m^f ar and m^ r according to phase space, the weight to be assigned to the configuration 

pear = l+ jfar = j- ig 



P + ^mu,m^r, mf{) = \ [f q P 2 (m lh m»r , mf) + fgP^mu, m^ , mf{ )j (2.2) 
while that for Z ncar = l~, Z far = l + is 



,near ™far\ , f r> /■„__ ^,ncar „_far\ 



1 



P- + (m u ,m?r,rn%n = \ [f q Pi(m n ,m^ ,mf) + f q P 2 (m u , m^ , mfi)\ . (2.3) 

In the former case, since the distinction between Z ncar and / far is lost in the data (except 
when interchanging them gives a point outside phase space), we must use eq. (2.1) with 
1+ = / ncar , l~ = Z far in the logarithmic factor of the KL-distance, i.e. the contribution is 

/ P + _(m // ,m° / car ,m^ ar |T) + P_+(m H , m% r , m°, car |T) \ 

l0 & T, 7 r^aT far! en , n 7 faT noar I a\ P +-\ m Ui m jl > m jl I 1 ) ■ \ 1A ) 

\P+- (mu , m*f ar , mf{ \ S) + P_+ (m« , m*f , m*f ar 1 5) y J J 

Similarly from the configuration P ear = l~ , Z far = l + we get the contribution 
/ P-+(mu, m?, car , mfrlT) + P+-(m lh m% r , m?, car |T)\ 

l^n- J' I / 1 i- V »' \> p , ^.near far^N / c\ 

i0g "5 7 T^aT fan ex | p 7 faT nearly P -+{ m ll,m j i ,m jt \1 ) . (2.5) 

y P- + (mu , m*jf ar , mjf 1 5) + P + _ (mjj , m*f , m°f ar 1 5) y 

Denoting the sum of these two contributions at the ith phase space point by KLj (T,S), 
and summing over M such points, we have as M — > oo 

M1 ° gfl AT, (2.6) 



E,KLi(T,S) 



which is the Monte Carlo equivalent of eq. (1.2). Results for R = 1000 and M = 5 x 10 7 
are shown in table 4. By comparing with tables 1-3, we see that, as might be expected, 
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the three-dimensional analysis achieves a discrimination that is better than that of a one- 
dimensional analysis applied to any single invariant mass distribution. This could be 
particularly useful in difficult cases like that of distinguishing between SFSF (SUSY) and 
FVFV (UED). 
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(a) 
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Table 1: The number of events needed to disfavour the column model with respect to the row 
model by a factor of 0.001, assuming the data to come from the row model, for the m\ distribution: 
(a) mass spectrum I and (b) mass spectrum II. 



(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 


SFSF FVFV FSFS FVFS FSFV SFVF 


SFSF 


oo 


1059 
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SFSF 


oo 


3006 958 
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761 
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FVFV 
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oo 
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FVFV 
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FSFS 


278 
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FSFS 
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FVFS 
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SFVF 


1141 
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4259 
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Table 2: As in 


table 1, 


for the 


m| i+ distribution. 








(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 


SFSF FVFV FSFS FVFS FSFV SFVF 
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Table 3: As in 


table 1, 


for the 


m|;_ distribution. 








(a) 


SFSF FVFV FSFS FVFS FSFV SFVF 


(b) 
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Table 4: As in table 1, for the combined three-dimensional distribution. 
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